Low-Energy Photodisintegration of the Deuteron and Big-Bang Nucleosynthesis 
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The photon analyzing power for the photodisintegration of the deuteron was measured for seven 
gamma-ray energies between 2.39 and 4.05 MeV using the linearly polarized gamma-ray beam of the 
High-Intensity Gamma-ray Source at the Duke Free-Electron Laser Laboratory. The data provide 
a stringent test of theoretical calculations for the inverse reaction, the neutron-proton radiative 
capture reaction at energies important for Big-Bang Nucleosynthesis. Our data are in excellent 
agreement with potential model and effective field theory calculations. Therefore, the uncertainty 
in the baryon density S7s/i 2 obtained from Big-Bang Nucleosynthesis can be reduced at least by 
20%. 

PACS numbers: 25.20.-x, 24.70.+S, 27.10.+h, 21.45. +v 
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Big-Bang Nucleosynthesis (BBN) is an observational 
cornerstone of the hot Big-Bang (BB) cosmology. Ac- 
cording to the neutron(n)-proton(p) capture reaction 
p(n,-f)d with a deuteron (d) and a 2.225 MeV 7 ray in 
the exit channel is of special interest, because the BB 
abundance of deuterium provides direct information on 
the baryon density in the early universe at times be- 
tween about 0.01 and 200 seconds after the BB. Know- 
ing accurately the n-p capture cross section in the en- 
ergy range from 25 to 200 keV in the center-of-mass 
(cm.) system and using the experimental value for the 
primeval deuterium number density (D/H) p would 
allow for an accurate determination of the baryon den- 
sity flgh 2 (h is the Hubble constant in units of 100 
km/s/Mpc). From fish 2 one can predict the abundances 
of the three light elements 3 He, 4 He, and 7 Li. Accord- 
ing to [l| , the 10% uncertainty in the deuterium- inferred 
baryon density Q B h 2 = 0.019 ± 0.002 comes in almost 
equal parts from the (D/H) measurements and theoreti- 
cal uncertainties in predicting the deuterium abundance. 
For the latter, the knowledge of the n-p capture cross 
section is of crucial importance. Unfortunately, there 
is a near-complete lack of data at energies relevant to 
BBN. Aside from thermal energies, data exist only at 
n-p cm. energies of 275 keV and above. As a conse- 
quence, the ENDF-B/VI evaluation has been used Q 
in the BBN energy range. This evaluation is normalized 
to the high-precision thermal n-p capture cross-section 
measurements. The 5% uncertainty that is assigned in 
this approach contributes a significant fraction to the un- 
certainty in the baryon density and consequently in the 
abundances of the light elements produced in BBN. 

Very recently, with the precision results from WMAP 
(Wilkinson Microwave Anisotropy Probe) for the Cos- 
mic Microwave Background (CMB) and its anisotropics 
an independent and even more accurate result became 



available: n B h 2 = 0.0224 ± 0.0009 HQ. The compari- 
son of the baryon density predictions from BBN and the 
CMB is a fundamental test of BB cosmology ■ Any de- 
viation points to either unknown systematics or the need 
for new physics. Therefore, it is of crucial importance to 
reduce the uncertainty in fi^/i 2 obtained from BBN. As 
stated above, 50% of the uncertainty is due to the un- 
certainty in the n-p capture cross section in the energy 
range of interest. 

Recently, effective field theory approaches 0, have 
provided accurate results for the time-reversed reaction 
7-d -> n-p from threshold (2.225 MeV) to about 10 MeV 
incident 7-ray energy. The work described in this Let- 
ter was motivated by these new theoretical results and 
also earlier nucleon-nucleon potential model based cal- 
culations |lC| in the 7-ray energy range important to 
BBN (E 1 '= 2.25-2.43 MeV). Aside from the 7-d cross 
section, and the n-p capture cross section inferred via 
"detailed balance," these calculations predict results for 
other observables as well which are related to the cross 
section, but are in principle experimentally easier to mea- 
sure with high accuracy than the n-p capture cross sec- 
tion itself. The aim of this work is to provide an alter- 
native method of determining the accuracy of theoretical 
models in predicting the n-p capture cross section in the 
energy range of interest for BBN. Potentially, this could 
lead to a considerably smaller uncertainty in Qsh 2 ob- 
tained from BBN. 

We measured the analyzing power S(90°) for the 
2 H(t*, np) reaction with linearly polarized 7 rays at 9 = 
90° (lab) for seven energies between E 1 = 2.39 and 
4.05 MeV. This energy range corresponds to n-p cm. 
energies of 165 keV to 1.83 MeV, i.e., the present experi- 
ment includes for the first time data in the upper energy 
range of interest to BBN. The analyzing power £(#) is 
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FIG. 1: Schematic of experimental setup. The 7-ray beam 
is perpendicular into the page and the 7-ray polarization is 
nominal in the <f> = plane. 



defined as 

U a(6»,0 = O°)+o-(0,0 = 9O°) / a + 6sin 2 6>'.f 
where the differential cross section a(9, </>) is given by 



bsm z 9 1 
■ x ■ 



a + 6 sin 2 9[1 + cos 20] 



Here, 6 is the polar angle, 4> is the azimuthal angle, and 
/ is the degree of linear polarization of the incident 7- 
ray beam. The quasi-monoenergetic and linearly polar- 
ized 7-ray beam was produced by Compton backscat- 
tering of relativistic electrons from 670 nm free-electron 
laser (FEL) photons at the High-Intensity Gamma-ray 
Source (HIGS) located at the Duke University Free- 
Electron Laser Laboratory. The electron energy in the 
electron storage ring was varied between E e = 300 and 
375 MeV to generate 7-ray beams of energy between 2.39 
and 4.05 MeV. At a distance of 75 m from the electron- 
FEL-photon collision point the collimated 7-ray beam 
of 2.6 cm diameter struck a 4 cm diameter and 6 cm 
long deuterated liquid scintillator (C6D12, Nuclear Enter- 
prises NE232) contained in a thin-walled glass container 
and viewed by a photomultiplier tube (PMT). The axis 
of the scintillator-PMT arrangement coincided with the 
axis of the incident 7-ray beam. The average 7-ray flux at 
the location of this deuterated scintillator target (DST) 
was 5 x 10 5 7/s. The 7-ray beam was monitored with a 
140% HPGe detector positioned downstream of the ex- 
perimental setup. Aside from low-energy 7-ray sources 
the "natural" 7-ray lines at E 1 — 1461 keV ( 40 K) and 
E 1 — 2614.5 keV ( 208 T1) served as convenient online cali- 
bration sources throughout the course of the experiment. 
The energy spread AE/E of the 7-ray beam varied be- 



tween 2.3% FWHM at E 1 = 2.39 MeV to 2.9% FWHM 
at E 1 = 4.05 MeV. 

The experimental setup is shown schematically in 
Fig. ^ Neutrons from the deuteron breakup reaction 
were detected by four Bicron 501A liquid scintillator de- 
tectors, 2" in diameter and 2" in length, viewed by a 
PMT. We used four detectors rather than two to increase 
the efficiency of our experimental setup. Two neutron 
detectors were mounted at 9i a b = 90° in the plane of 
the 7-ray polarization (nominally the horizontal plane) 
on opposite sides of the incident 7-ray beam (cf> = 0° 
and 180°). The other two detectors were mounted at 
Slab = 90° in the perpendicular plane (</> = 90° and 
270°). The center-to-center distance between the DST 
and the neutron detectors was 17 cm. The protons from 
the deuteron breakup in the DST gave the start sig- 
nal for a neutron time-of-ffight measurement between the 
DST and the neutron detectors. Neutron-gamma pulse- 
shape discrimination (PSD) techniques were applied to 
distinguish the events of interest from the overwhelming 
background produced in the neutron detectors by Comp- 
ton scattering from the DST. Two-dimensional spectra 
of pulse height in the DST versus neutron time-of-ffight 
were created for the four neutron detectors used in the 
present experiment. Time-of-ffight and proton recoil en- 
ergy spectra are shown in Figs. [21 and [3] for an incident 
gamma-ray energy of 4.05 MeV. Our experimental tech- 
niques cannot be extended to much lower 7-ray energies 
than already achieved in the present experiment. At E 7 
= 2.39 MeV both the proton and neutron energies were 
only 90 keV compared to the more comfortable value of 
about 900 keV at E 1 = 4.05 MeV. Liquid scintillator de- 
tectors are not commonly employed to detect neutrons 
and protons at energies of less than 500 keV. However, 
other types of detectors lack the fast timing character- 
istics and efficiencies that are crucial for obtaining 7-d 
data in the energy range between E 1 = 2.4 and 3 MeV. 

In order to cancel instrumental asymmetries in our ex- 
perimental setup, we either rotated the neutron detec- 
tors which were mounted on a ring centered and posi- 
tioned perpendicular to the 7-ray beam axis through 90° 
(counter clockwise), or we interchanged the detectors, 
i.e., the detectors in the horizontal plane were moved to 
the vertical plane and vice versa. Within statistical un- 
certainties either procedure gave consistent results for the 
asymmetry e, which was calculated from the formula e = 
[a — l)/(a+l). For the rotation procedure we defined a 
asai_ 2 = [(JVf R N^ L )/{N^ U NY U )] 1/2 for detector pair 
1-2, and as a 3 - 4 = [{N^ L N^ R ) / N^ N^ )] 1 ' 2 for de- 
tector pair 3-4. Similarly, for the interchange procedure 
we have «i_ 2 = [{N^ r N^ r )/{NY u NY U )} 1/2 for detec- 
tor pair 1-2, and a 3 _ 4 = [{N^ L N^ L )/{NX D N^ )] 1 ' 2 
for detector pair 3-4. Here, Nf IR (Nl 1L ) refer to the neu- 
tron yields detected with detector i positioned in the hor- 
izontal plane to the right (left) side of the incident 7-ray 
beam, and N^ 17 (Ny D ) refer to the neutron yields de- 
tected with detector i positioned in the vertical plane in 
the up (down) position. 
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FIG. 2: Neutron time-of-flight spectrum between the deuter- 
ated scintillator and a neutron detector for the reaction 7- 
d — > 71-p at _E 7 = 4.05 MeV. Time increases from left to right. 
The dominant peak is due to the neutrons of interest. The 
small peak is due to 7 rays leaking through the PSD cut. This 
leakage is smaller than 0.1%. 
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FIG. 3: Proton recoil energy spectrum in the deuterated scin- 
tillator (DST) at £ 7 = 4.05 MeV. The small peak near chan- 
nel 450 is due to electrons generated via Compton scattering 
to the neutron detector. 



Based on the geometry of the undulator magnets used 
to produce the FEL photons, the photon polarization 
should be linear and of magnitude 1.0. Furthermore, 
the polarization plane should coincide with the horizon- 
tal plane in the laboratory. Using the polarization de- 
pendent formulas for inverse Compton scattering this 
should result in a linear 7-ray polarization in the horizon- 
tal plane of / = 1.0 in the photon and electron energy 
range of interest for the present experiment. However, 
the optical cavity mirrors used to produce the FEL pho- 
tons of 670 nm were optically active, causing a rotation 
of the polarization plane of the FEL photons and con- 



sequently of the resulting 7-ray beam. The outcoupled 
FEL light was used to verify the linear polarization of 
1.0 and to determine the tilt angle of the polarization 
plane. However, there is no guarantee that the tilt an- 
gle of the polarization plane of the outcoupled light (i.e., 
the transmitted light through the mirror opposite to the 
7-ray beam direction) is in perfect agreement with the 
tilt angle of the FEL photon polarization inside of the 
optical cavity where the electron-photon collision takes 
place. Therefore, the tilt angle of the 7-ray polarization 
plane relative to the nominal horizontal plane was deter- 
mined from the measured asymmetry e of the Compton 
scattered 7 rays by setting the PSD gate on the 7 rays 
in the neutron detectors and by selecting the appropri- 
ate pulse height gate (due to 7-ray scattering from elec- 
trons through 90°) in the DST. This asymmetry e was 
determined simultaneously with the one for the breakup 
neutrons from the j-d reaction. In order to extract the 
tilt angle from the measured 7-ray asymmetry data, the 
effective analyzing power of our apparatus for Compton 
scattering from electrons was calculated via Monte-Carlo 
simulation using the Klein-Nishina formula. The polar- 
ization in multiple Compton scattering was treated ex- 
actly. The average tilt angle of the 7-ray polarization 
plane was found to be (13.7 ± 0.2)° in upward direction 
relative to the horizontal laboratory plane. This value 
is about 2° larger than the polarization tilt angle of the 
outcoupled FEL photons. 

The neutron asymmetry data from the j-d — > n-p re- 
action were corrected for finite geometry and multiple- 
scattering effects via extensive Monte-Carlo simulations 
of the experimental setup, using the tilt angle of the 7-ray 
polarization determined above and the j-d cross section 
and analyzing power calculations of Arenhovel 01 which 
are based on the Bonn nucleon-nucleon potential model 
|ll|. The use of an active deuterium target makes our 
data practically insensitive to multiple 7-ray scattering 
(i.e., Compton scattering off electrons) in the DTS before 
the 7-d — > n-p reaction of interest is taking place. The 
light output produced by the recoil electrons generated 
in the Compton scattering process is considerably larger 
than the light output produced by the protons from the 
7-d — > n-p reaction. Therefore, multiple 7-ray scattering 
can be eliminated efficiently by setting a tight gate on 
the proton pulse height of interest. In contrast, multiple 
scattering of the neutrons from the 7-d reaction in the 
DST has to be taken seriously. Especially at the lowest 
7-ray energies employed in the present experiment, our 
constraint on the proton pulse height in the DST and 
our cut on the neutron time-of-flight did not eliminate 
multiple scattering events completely due to limitations 
of the detectors' energy and time resolutions. 

The data for the excitation function of the analyzing 
power S(90° lab) are shown in Fig.0|and listed in TableQ] 
The error bars include statistical and systematic uncer- 
tainties added in quadrature. At the higher energies the 
analyzing power is close to 1, i.e., the neutrons are emit- 
ted almost completely in the plane of the 7-ray polariza- 
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tion (electric dipole radiation El). At energies below E 1 
— 3 MeV, £(90°) decreases rapidly, i.e., the probability 
for neutrons to be emitted in the vertical plane (mag- 
netic dipole radiation Ml) increases with decreasing 7- 
ray energy. The curve shown in Fig. ^is the prediction of 
Arenhovel ^3] using the coordinate-space version of the 
Bonn nucleon-nucleon potential model The calcu- 



lation includes meson-exchange, isobar, and relativistic 
effects. Clearly, the model calculation is in very good 
agreement with the experimental data. Table [i] shows 
that the effective field theory approach of Chen and Sav- 
age |8j gives basically the same results as the potential 
model calculation of Arenhovel. 



TABLE I: Measured photon analyzing power E at 9i a t = 90° 
in comparison to theoretical predictions. 

E-y (MeV) 6cm. (deg) (keV) £ E Ar enhoevel Schen&S* 



2.39 


95.6 


165 


0.419 ±0.021 


0.461 


0.464 


2.48 


94.6 


255 


0.649 ± 0.019 


0.624 


0.631 


2.60 


94.0 


375 


0.745 ± 0.022 


0.760 


0.757 


3.02 


93.2 


795 


0.911 ±0.014 


0.902 


0.901 


3.22 


93.0 


995 


0.928 ±0.012 


0.925 


0.923 


3.52 


92.9 


1295 


0.953 ±0.012 


0.944 


0.942 


4.05 


92.8 


1825 


0.975 ±0.013 


0.959 


0.958 



TABLE II: Ml (s-wave) contribution S to the total 7-d cross 
section obtained from the present £(#) data in comparison to 
the predictions of Arenhovel, Chen and Savage, and Rupak. 



(MeV) 


S 

this experiment Arenhovel Chen & Savaj 


;e Rupak 


2.39 


0.675 ±0.019 


0.662 


0.622 


0.627 


2.48 


0.448 ± 0.020 


0.468 


0.459 


0.458 


2.60 


0.339 ± 0.026 


0.328 


0.320 


0.317 


3.02 


0.128 ±0.019 


0.141 


0.139 


0.135 


3.22 


0.104 ±0.017 


0.108 


0.109 


0.104 


3.52 


0.069 ±0.017 


0.080 


0.083 


0.079 


4.05 


0.037 ±0.019 


0.057 


0.061 


0.056 



As shown in detail by Schreiber et al. [l2j the 7-d ana- 
lyzing power data £(6>) at low energies can be used to de- 
termine the relative Ml and El strengths of the 7-d cross 
section. Based on the present S (90° lab) data, Table UTI 
gives the calculated Ml contribution to the 7-d cross sec- 
tion in comparison to the effective-field theory calcula- 
tions of Chen and Savage, Rupak 0, and the nucleon- 
nucleon potential-model calculation of Arenhovel. In the 
energy range most important for BBN our results are in 
excellent agreement with the theoretical predictions, es- 
pecially with the calculations of Arenhovel. Fig. [S] shows 
the calculated total 7-d cross section of Chen and Sav- 
age as well as the associated Ml and El contributions 
in comparison to the Ml contribution determined in the 
present work (dots) and in the earlier work of Schreiber 
et al. at E 1 = 3.58 MeV (triangle). 

In summary, the first experimental test of theoretical 
models used to calculate the n-p capture cross section in 
the energy range of importance to BBN reveals almost 
perfect agreement with experimental information derived 
from analyzing power data for the reverse reaction 7- 
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FIG. 4: Excitation function of the photon analyzing power £ 
for the reaction 7-d — > n-p at 0i a b = 90° in comparison to the 
theoretical prediction of Arenhovel. 



d — > n-p. This observation lends substantial credibility 
to the theoretical models also in the presently not tested 
7-ray energy range from 2.25 to 2.38 MeV, i.e., for n-p 
cm. energies between 25 and 155 keV. Work is planned 
to reduce the uncertainty of the described measurements 
from its present 3% uncertainty at 2.39 MeV in deter- 
mining the dominant Ml contribution to the 7-d cross 
section to 1.5% and to extend the measurements to n-p 
cm. energies as low as 25 keV. 

We conclude that the ±5% uncertainty used in Q is 
a very conservative estimate for the uncertainty of mod- 
ern theoretical approaches available for calculating the 
n-p capture cross section in the energy range relevant to 
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FIG. 5: Data for the Ml contribution to the ■y-d total cross 
section in comparison to the theoretical prediction of Chen 
and Savage (dashed curve). The dashed-dotted curve repre- 
sents the El contribution and the solid curve in the total 7-d 
cross section from Chen and Savage. Note the logarithmic 
scales for both a and E y . 



BBN. The uncertainty quoted in Ref. pj can be reduced 
by at least 20%. The planned improvements of our mea- 
surements are expected to provide an even more accurate 
test of the calculated n-p capture cross section. There- 
fore, this cross section will play a small role in the overall 
uncertainty of the baryon density in the early universe 
as determined in Ref. The BBN approach compares 
favorably with the very recent CMB based method of 
determining n B h 2 from the WMAP data [EE!- 
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